Network topology and collapse of collective stable chaos 
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Collective stable chaos consists of the persistence of disordered patterns in dynamical spatiotem- 
poral systems possessing a negative maximum Lyapunov exponent. We analyze the role of the 
topology of connectivity on the emergence and collapse of collective stable chaos in systems of cou- 
pled maps defined on a small- world networks. As local dynamics we employ a map that exhibits a 
period-three superstable orbit. The network is characterized by a rewiring probability p. We find 
that collective chaos is inhibited on some ranges of values of the probability p; instead, in these re- 
gions the system reaches a synchronized state equal to the period-three orbit of the local dynamics. 
Our results show that the presence of long-range interactions can induce the collapse of collective 
stable chaos in spatiotemporal systems. 
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I. INTRODUCTION 

The study of dynamical processes on nonuniform or 
complex networks has attracted much attention (Motter 
ct. al., 2006; Newman, Barabasi and Watts, 2006). In 
particular, the influence of the network topology have 
been investigated on phenomena such as synchroniza- 
tion (Barahona and Pecora, 2002), pattern formation 
(Cosenza and Tucci, 2001), spatiotemporal intcrmittency 
(Cosenza and Tucci, 2002), phase-ordering and critical 
phenomena (Dorogovtsev, Goltsev and Mendes, 2008), 
among others. 

In this context, there has been interest in the research 
of the phenomenon of nontrivial collective behavior in 
dynamical networks. This phenomenon arises when dis- 
similar forms of temporal evolution of macroscopic and 
microscopic quantities coexist in a system (Manrubia, 
Mikhailov and Zanette, 2004). Two different manifes- 
tations of nontrivial collective behavior have been inves- 
tigated. On one hand, there is emergence of order in the 
time evolution of the dynamics of macroscopic quanti- 
ties of a system of chaotic elements. For example, the 
spatial average of the states of the elements may be pe- 
riodic in time, while the evolution of the individual dy- 
namics of these elements is chaotic and desynchronized. 
This kind of collective behavior has been widely studied 
(Kaneko, 1990a; Kaneko, 1994; Chate and Manneville, 
1992; Cosenza and Gonzalez, 1998). The counterpart of 
this phenomenon has also been observed: the persistence 
of collective chaos or temporal disorder at the macro- 
scopic level in systems of interacting elements whose in- 
dividual dynamics is periodic. This nontrivial behav- 
ior has been called collective chaos (Politi, Livi, Oppo 
and Kapral, 1993; Kapral, Livi, Oppo and Politi, 1994; 
Kapral, Livi, Oppo and Politi, 1997) and is one of the 
least understood emergent phenomena in complex sys- 
tems. 

Stable collective chaos consists of an irregular behav- 
ior that can not be described by the presence of repellers 
in the phase space of the system, producing, as a conse- 



quence, the divergence of nearby trajectories. Moreover, 
in such systems, the duration time of the transient regime 
can scale exponentially with respect to the size system, 
and the final stable attractor can not be reached in the 
practical sense for large enough systems (Crutchfield and 
Kaneko, 1988; Kaneko, 1990b; Politi et al, 1993; Cec- 
coni, Livi and Politi, 1998; Bagnoli and Cecconi, 2001). 

This type of irregular collective behavior that occurs in 
systems of coupled periodic dynamic elements is different 
from the phenomenon of transient chaos, which is a truly 
chaotic regime with a finite lifetime, and is characterized 
by the coexistence of stable attractors and non-attractive 
chaotic sets (repellers) in the phase space of a system (Tel 
and Lai, 2008). In such systems, a generic initial config- 
uration typically produces an irregular trajectory until it 
suddenly collapses into a non-chaotic attractor (Wacker- 
bauer and Showalter, 2003). In the case of stable col- 
lective chaos, the transient collective behavior, which is 
usually considered irrelevant, is chaotic and statistically 
stationary. The behavior of the system in the transient 
regime can not be distinguished from a typical chaotic 
behavior. For this reason, we speak of the existence of 
supcrtransicnts, since the chaotic collective behavior is, 
even for a system of moderately small size, the only be- 
havior observable in practice. This phenomenon was first 
reported in a network of coupled chaotic maps (Crutch- 
field and Kaneko, 1988; Kaneko, 1990b). This system 
may exhibit, for a system of 128 elements and taking 
into account the speed and precision of the computer, a 
characteristic time of the supertransitory regime of about 
10 64 years (Crutchfield and Kaneko, 1988). This result 
has a deep impact on some physical processes that are not 
yet fully understood, such as turbulence in fluids or the 
observation of aperiodic behaviors in complex systems. 
Such processes could correspond, from a strictly mathe- 
matical point of view, to a transient state. In practice, 
we may never reach the final regular behavior. Thus, the 
truly stationary, observable behavior will be the super- 
transitory regime of collective chaos (Bagnoli and Cec- 
coni, 2001). 



2 



Transient spatiotemporal chaos has been studied in 
reaction-diffusion continuous systems, such as Gray- 
Scott's equations (Wackcrbaucr and Showaiter, 2003; 
Wackerbauer and Kobayashi, 2007; Wackerbauer, 2007; 
Yonker and Wackerbauer, 2006), where it has been found 
that the spatial boundary conditions can induce the col- 
lapse of transient chaos towards a fixed point. On the 
other hand, long transients appear in networks of model 
neurons when the number of connctions per neuron is 
small (Zillmer, Brunei and D. Hansel, 2009). 

In this article we investigate the influence of the topol- 
ogy of the connectivity in the process of collapse of col- 
lective chaos in networks of interacting dynamical ele- 
ments. Recently, intensive research has been performed 
on the theory and applications of small-world networks 
(Newman, Barabasi and Watts, 2006). They constitute 
a class of networks with a high degree of local cluster- 
ing and a small characteristic length between any two 
elements (Watts and Strogatz, 1998). It has been shown 
that small- world networks describe many natural and ar- 
tificial networks. By varying a parameter, small-world 
networks can be continuously tuned between ordered, de- 
terministic lattices and completely random networks. In 
this paper we consider a coupled map model defined on a 
small- world network that exhibits the phenomenon of col- 
lective chaos and show its collapse as a result of the vari- 
ation in the connection topology of the system. Because 
of their discrete spatial nature, coupled map systems are 
especially appropriate for investigating dynamical pro- 
cesses occurring in heterogeneous media or in complex 
networks. 



II. THE MODEL. 

We use the algorithm of construction of small-world 
networks originally proposed by Watts and Strogatz 
(Watts and Strogatz, 1998). We start from a ring with 
N nodes, where each node is connected to its k nearest 
neighbors, k being an even number. Then each connec- 
tion is rewired at random with probability p to any other 
node in the network, to avoid self-connections. After the 
rewiring process, the number of elements coupled to each 
site - which we call neighbors of that site - may vary, but 
the total number of links in the network is constant and 
equal to Nk/2. It is assumed that all links are bidirec- 
tional. The condition log TV <C k <C N is employed in 
the algorithm to ensure that no node is isolated after the 
rewiring process, which results in a connected graph. For 
p = 0, the network is completely regular, while for p = 1, 
the resulting network is completely random. With this 
algorithm, a small- world network is formed for values of 
the probability in the range p £ (0.001, 0.2). 

The state of each node in the network can be charac- 
terized by a continuous variable, which evolves according 
to a deterministic rule depending on its own state and 
the states of neighboring nodes. We define a diffusively 



coupled map system on this network as 

4u = + [f{xi)-f{x\) 



(i) 



where x\ is the state of i-th map for the discrete time 
t, with i = 1,2, ... ,N; k l is the number of elements 
connected to the element i; v l is the set of neighbors of i; 
f(x) is a function that describes the local dynamics; and 
e is the coupling parameter. 

A map with minimal ingredients that presents collec- 
tive chaos is (Politi et al., 1993) 
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(2) 



where x G [0, 1], and the parameters of the map a, 6, and 
c are chosen such that the local dynamics converge to- 
wards a period three superstable orbit, corresponding to 
the points x\ = a — > X2 = ab — > x% = ab 2 . A superstable 
orbit has a Lyapunov exponent, A, with value A — > — oo. 
Here we shall set a = 0.1, b = 2.5, and c = 0. 

Note that the stability of the period-three orbit of the 
local map Eq. ([2]) implies the stability of the system given 
Eq. ([1]), whose maximum Lyapunov exponent turns out 
be negative for any value of e. Thus, the long-term evo- 
lution of the coupled system is confined to this periodic 
attractor. 

The collective behavior of the network can be charac- 
terized by the instantaneous mean field of the system, 
defined as 
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Similarly, the state of synchronization of the elements 
in the system Eq. ([T} can be characterized by the asymp- 
totic time average (a) of the instantaneous standard de- 
viations a t of the distribution of variables x\, defined as 
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A synchronized state corresponds to (a) = 0. In our 
case, the system Eq. (JTJ) can be synchronized in the 
period-three superstable orbit of the local map for certain 
values of parameters. 

In all calculations shown below, we use networks with 
size N = 10 4 , and assign random initial conditions x\ € 
[0, 1] in the system Eq. (JT]). 

Figure [T] shows the bifurcation diagram of the mean 
field, H t , and the average standard deviation (a) as func- 
tions of coupling parameter e for a network with k = 2 
and p = 0, which corresponds to a one-dimensional reg- 
ular lattice. 

In Figurc[T]we see that the macroscopic variable H t dis- 
plays chaotic behavior and chaotic bands in a wide range 
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FIG. 1: Bifurcation diagrams of the mean field Ht (top) and 
the quantity (a) (bottom) as functions of e, for the system 
Eq. (JTJ with k = 2 and p = (a one-dimensional lattice), 
size N = 10 4 . For each value of e, 10 3 values of Ht are 
plotted, after discarding 10 4 iterations. The quantity (a) was 
calculated as the average of 10 3 values of a t , after discarding 
10 4 iterations. 



of the parameter e, in spite of the existence of period- 
three superstable orbit in the local maps. The global 
attractor of the system must be a synchronized orbit 
equal to the superstable period-three of the local maps. 
However, we observe that the quantity (a) 0, indicat- 
ing that the system does not reach a synchronized state 
in its long-term evolution. The convergence time to get 
the synchronized periodic asymptotic state is extremely 
large. Figure [1] shows evidence of the phenomenon of 
collective stable chaos occurring in this system. 

To study the effects of long-range interactions in the 
phenomenon of stable chaos in the system Eq. (|T|), we 
calculate the quantity (er) as a function of the rewiring 
probability p, as shown in Figure [2] There is a criti- 
cal value of the probability p above which the system 
reaches a synchronized state, characterized by (a) — > 0. 
This synchronized state corresponds to the superstable 
orbit of superstable period-three of the local map Eq. © . 
Thus, the presence of long-range interactions associated 
to an increase of p contributes to the collapse of collective 
stable chaos and to the emergence of synchronization in 
this system. 

The transient nature of collective chaos observed in 
the system Eq. ([T]) implies that the chaotic behavior ob- 
served for certain ranges of p in Fig. [5] will eventually 
converge towards the stable period-three attractor of the 
local map. However, the regime of collective chaos is, in- 
deed, supertransient for those parameter values. We have 
verified that, even for a small system consisting of N = 1 
maps, it requires about 10 10 iterations to reach the syn- 
chronized period-three state. Thus, in practice, collective 
stable chaos becomes the actual observable macroscopic 
state for those values of parameters. 
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FIG. 2: The quantity (a) as a function of the reconnection 
probability p, for fixed e = 0.384, k = 10, and N = 10 4 . The 
quantity (a) is the average of 10 2 values of at, after discarding 
3 x 10 4 iterations. Error bars correspond to the standard devi- 
ations of 10 realizations of initial conditions for each value of 
p. The numbers indicate the regions of the parameter p cor- 
responding to (1) regular network, (2) small-world network, 
and (3) random network. 



III. CONCLUSIONS 

The existence of a supertransitory regime such as col- 
lective stable chaos can not be distinguished from a typ- 
ical stable chaotic behavior in a system. For this rea- 
son the supertransitory states are, even for a moderately 
small system, the only real observable behavior. 

Figure [5] shows that the system can reach the synchro- 
nized state of period three by changing the topology of 
the network. The presence of long-range interactions be- 
tween the nodes, favored by increasing the rewiring prob- 
ability p, facilitates the transmission of information and 
the possibility of synchronization in small- world and ran- 
dom networks, and therefore contributes to the collapse 
of collective stable chaos. 

Our results indicate that the persistence of collective 
chaos in spatiotcmporal systems depends on topological 
properties of the spatial substrate and, therefore, are not 
robust properties. The results also suggest that non- 
trivial collective behavior observed in many complex sys- 
tems can be modulated or controlled through changes in 
their geometric properties. 
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